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Abstract 

A multiple-input single-output (MISO) wiretap channel model is considered, that includes a multi- 
antenna transmitter, a single-antenna legitimate receiver and a single-antenna eavesdropper. For the 
scenario in which spherical uncertainty for both the legitimate and the eavesdropper channels is included, 
the problem of finding the optimal input covariance that maximizes the worst-case secrecy rate subject 
to a power constraint, is considered, and an explicit expression for the maximum worst-case secrecy rate 
is provided. 
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I. Introduction 

Wireless physical (PHY) layer based security approaches exploit the physical characteristics of the 
wireless channel to enhance the security of communication systems. The most basic physical layer model 
that captures the problem of communication security was proposed in m. Later, the Gaussian scalar 
wiretap channel was studied in ||2l. Recently, the secrecy capacity of multi-antanna wiretap channels has 
been studied in 131, H. The above work assumes that perfect channel state information (CSI) on the 
legitimate and the eavesdropper channels. When ellipsoidal channel uncertainty is included, the worst- 
case secrecy rate of multiple-input single-output (MISO) cognitive radio network was studied in ||5l. In 
151 . the worst-case secrecy rate maximization problem is converted to a quasi-convex problem by using 
the S-procedure to express the channel uncertainty constraint in a linear matrix inequality (LMI) form. 
In l6l, a MISO channel in the presence of multiple eavesdroppers, each equipped multiple antennas, was 
considered. Perfect CSI as well as channel uncertainty were considered in 0. As in 121, the problem 
was converted to a semidefinite programming (SDP) problem by using the S-procedure to express the 
channel uncertainty constraint in LMI form. 

In this paper, we consider the same MISO wiretap scenario as in Q, 0, except that we consider a 
special case of the channel uncertainty considered in Q, l6l, i.e., spherical uncertainly. For this case, we 
derive an explicit expression of optimal input covariance that maximizes worst-case secrecy rate subject 
to a power constraint. In particular, the solution is obtained via finding the eigenvalues of a known 6-by-6 
matrix. The advantage of such an explicit solution lies in: 

1) Independent of the number of antennas the problem leads to a 6-by-6 matrix whose entries are 
obtained by direct and simple calculation from the estimated channel values (see Eq. ([8]l and ([281)). 
The computation time for finding the eigenvalues of such 6-by-6 matrix is not affected by the size 
of problem. In contrast, for the existing methods (Q, 0), the computation time for the iteration 
algorithm is affected significantly by the size of problem. 

2) The structure of the optimal input covariance matrix is explicitly given (see Eq. (fT2l)). 

3) A (simple) necessary and sufficient condition for a positive worst-case secrecy rate is explicitly given 
(see Theorem [1]). 

Notation- Upper case and lower case bold symbols denote matrices and vectors, respectively. Super- 
scripts *, T and f denote respectively conjugate, transposition and conjugate transposition. A ^ means 
that the matrix A is Hermitian positive semi-definite. \a\ denotes the absolute value of a, while ||a|| 
denotes Euclidean norm of the vector a. eig(A) denotes the eigenvalues of the matrix A. 
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II. System Model and Problem Statement 

Consider a Gaussian MISO wiretap channel, which includes a transmitter is equipped with nj- antennas, 
and a legitimate receiver and an eavesdropper, each equipped with a single antenna. The received signals 
at the legitimate receiver and the eavesdropper are respectively given by 

UR = hjjx + VR, and ue = h^x + ve (1) 

where x is the x 1 transmitted signal vector with zero mean and covariance matrix Q = E{xx^}; 
h.R, h^; are respectively channel vectors between the transmitter and legitimate receiver, and between the 
transmitter and eavesdropper; vr, ve are the noises at the legitimate receiver and the eavesdropper with 
zero means and unit variances, respectively. 

We consider the scenario in which spherical channel uncertainty, i.e., Q §18.2.4] 

hR£nR = {hR\\\hR-hR\\<eR}, (2) 

hE enE = {hE\\\hE-hE\\<eE} (3) 

where h.R and h^; are the estimated values, eR and e^; are the estimated error bounds. Let P be the 
transmitted power budget. We investigate the maximum worst-case secrecy rate defined as |I3 

Rs = max min i?s(Q, h^, h^;) (4) 

Q^O,Tr(Q)<P hflenn.hEGHE 

where Rs{Q,h.R,h.E) = log(l+h^Qh/i;)— log(l+h^Qh£;). We aim to find the optimal input covariance 
matrix Q*. A necessary condition to ensure Rs > is ^ ^r, i.e., 

Ilhfill > eR. (5) 

We assume that the condition dSjl holds in this paper. 

Since log(-) is an increasing function, (0]) is equivalent to 

l + h^RQhR 

max mm . (6) 

QbO,Tr(Q)<P h„,eOfl,hE6nB 1 + h^Qh^ 

Let T* be its optimal objective value. Certainly, r* > 1 is equivalent to Rg > 0. 

III. Explicit Solution of Optimal Input Covariance Matrix 

One can show that if r* > 1 {Rs > 0), then there exists a rank one matrix, to be denoted by Q*. The 
proof can be found in [61]. With this, we assume r* > 1 and find the rank one Q*. If such rank one Q* 
achieves Rs > 0, Q* is indeed the solution of the problem. Otherwise, Rg = 0. 
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Since r* > 1, it is easy to verify that Tr(Q*) = P. Then, we let Q = Puu^" and rewrite Q as 

1 + P min |u^^h/j|^ 

1 I D 1 12 - 

||u||2=i 1 + -P max |uTh£;p 



Let us denote 



Ihiill ||h£;|| 

- II II , and zq = max(c, \/l — ) . (8) 
||n£;|| 



Let z* be the solution of 

l + h{z- cf 



max 

20<2<1 



1 + a(rz — \/l — r^\/l — + d) 
The main result of this paper is given in the following theorem. 

Theorem 1: 

i) The sufficient and necessary condition for Rs > is that 



Vl -r2 Z0+ {r- ^/bja ) ^/l - z'^ > 
\/l-r'^^/l-z^ - {r-^Jhfa )2o > c^/hja + d 



(9) 



(10) 



Vl - r2 zo + (r - ) ^1 - 4 < 
or < y — . (11) 

w 1 - + {^/h/a- r)2 > c^Jhla + d 

ii) A55Mme f/zaf f liOD or f liil ) holds. In other words, t* > 1. T/ie solution of (O 



can Z^e seen from f li2D f/iaf u* a linear combination ofhji and Iie- The maximum worst-case 
secrecy rate is 

1 + biz* - cf 

Rs = ^og ^ . — 9- (13) 

1 + a{rz* - Vl - rVl - z*^ + fi) 

Remarks: Let us gain some insight into the conditions for Rs > 0. Obviously, (fTOl) and ([TT]) are both 
independent of P. In other words, if neither of these two conditions holds, Rs> Q cannot be achieved even 
when infinite power is used. Roughly speaking, this may occur if the estimation is not good enough (c, d 
is relatively large), or the legitimate channel is worse than the eavesdropper channel, i.e., ||hi?|| < \\h.E\\- 

Proof of Theorem [l] 



August 2, 2011 



DRAFT 



4 



For the optimal u, consider the two subproblems in ([7]): 

min |u^h/jp, (14) 

and max |u^h^|^. (15) 

Their optimal objective values can be obtained in closed form, given by ( | hji | — e/j ) ^ and (|u1^h£;|+e£;)^ 
respectively. Moreover, since r* > 1, it holds that |u^h/j| > cr. More details can be found in Appendix 
lAl From this result, ([7]) is equivalent to 

l + P(|uthj?|-eH)^ 
max_ - — I — . (16) 

||u|p = l, |uthH|>en 

We can reduce ([T6l ) to a problem of one variable. To this end, let z = |u1'h^|/||h/j ||. The domain of z is 

< ^ < 1, (17) 

W^rW 

which follows from Cauchy's inequality and the constraint |uth/j| > cr. It can be seen that the condition 
^ ensures a non-empty domain of z. Obviously, for fixed z (i.e., |u1'h/{| is fixed at ||h/j||z), ju^'h^;! 
should be minimized. Keeping this in mind, let 

=min |u^h£;| (18) 

u 

s.t. ||u|p = 1, and |u'''hj:j| = ||hR||z. 

Then, we rewrite (fT6ll as 

max — ^, , , , -75- s.t. ,,- ,, < z < 1. (19) 

z i + P{^{z) + €eY llhfill 

The function il^{z) is obtained in closed form using the following lemma. 

Lemma 1: Let a and b be (known) unit-norm vectors with r = |b^a| < 1, and < (7 < 1. The 
problem 

mm u' bbtu (20) 

u 

s.t. u^aa^u = g, and ||u|p = 1 

has an optimal objective value 



{r^/q — \/l — r^-^l — q)"^ if q > 1 — r 



2 

(21) 

if < 1 - 



Moreover, if q > 1 — r^, then the optimal u is given by 
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The proof is given in Appendix iBl The proof is similar to the one of Lemma 2 in |[8l. 

We assume that Iir and h^; are linearly independent (otherwise, the problem is reduced and simpler). 
Let a = h/j/||h/j||, b = h£;/||h£;||, r = |h^h/j|/(||h/i;|| ||hs||), q = z^. It follows from Lemma [U that 

r ||h£;||(rz — \/l — r^Vl — -z^ ) z> \/l — r'^ 
i{j{z) = I . (23) 

otherwise 
Note that in ( fT9l ). (|h/{|2: — eji)'^ is an increasing function of z. On the other hand, ip{z') = for any 
z' < \J\ — r"^ . Thus, it holds that z* > \J\ — r"^ (otherwise, there exists < z' < \/l — such that z' 
achieves a larger objective value. But this contradicts the optimality of z*). With this, we rewrite (O as 
do]). Here we replace the constraint 2 > c by z > c for convenience. Furthermore, from (l22l) . we obtain 

Gal). 

From the objective in the sufficient and necessary condition for i?^ > is that there exists a 
zq< z <\ such that \fh (z — c) > ^/a {rz — Vl — r^Vl — z'^ + d) , or equivalently 

max = \J \ — r'^\J \ — z"^ — {r — \fbfa )z] 

>c^/h/a + d. (24) 

Consider the optimization problem in the left hand side of ((24l) . The derivative g'i{z) = —z^/l — r^/ \/l — — 
(r — \JhJa ) is a strictly decreasing function of z. On the other hand, g'^iz) — )• —00 as 2; — 1. Thus, if 
9'\{zq) < 0, then zq is the solution; if g'i{zo) > 0, then the optimal z is the unique point with zq < z < 1 
and g'l (z) = 0. From this fact, we obtain (ITOl ) or (fTTl ). This completes the proof. 



A. Determining z* 

Note that the derivative g'{z) — )• —00 as z — )• 1. Thus z* is either zq or one of the feasible points (if 
any) with g'{z) = 0. We introduce the bijective transform 

2x 



1 + X2' 



< X < 1 



to rewrite Q as 



max 

X 



Fix) 



POX'^ + piX^ + P2X'^ + piX + po 



qox'^ + g'lX'^ + q2X^ + ^32; + ^4 



(25) 



(26) 



s.t. 



1 



Zq 



< x<l 



where po = l + bc^, pi = -Abe, p2 = 46 + 260^ + 2, qo = l + a(\/l - r^ + d)^, = 4ar(\/l - + 
(?2 = 2 - 2a + ear^ + 2a(f, qs = Aar{-Vl - r"^ + d), 94 = 1 + a{-y/l - + The solution x* 
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should be (1 — a/1 — Zq )/zq or one of the feasible points (if any) with F'{x) = 0. One can rewrite 
F'{x) = as an equation of six degree 

a^x^ + aix^ + 02^^ + a^x^ + a^x"^ + a^x + og = (27) 

where oq = piQ'o - Po^i, oi = 2p2Q'o - 2pog2, 02 = Spi^o + PiQi - Piq2 - ^Poqs, ^3 = 4pogo + 2pigi - 
2pig3 - 4g4Po, 04 = Spoi?! + Pig'2 - P2q3 - 3qipi, 05 = 2pog2 - 2g4P2, ^6 = Po93 " 94^1- AH the six 
roots of the equation (l27l) equal the eigenvalues Xi of the following 6x6 matrix fTOl Ch. 6] 



G 



ao ao flo f'-o (^o f^o 



(28) 



eisfG 



(29) 



I5 Osxi 

Thus, the optimal x satisfies (if there exists several feasible candidates for x*, the best one is chosen) 

X* G {■ 
Once X* is obtained, we have 

= T^2- (30) 

1 + x*^ 

IV. Numerical Simulations 

In this section we provide some examples to illustrate the result. We focus on how to obtain the explicit 
solution. For more numerical simulations on different estimated channel values, e.g., (e/?, es) please refer 
to El and in. 

First consider a MISO wiretap channel with = 4 antennas. We set {eji,eE) = (10~^,10~^), 
P = 5 dB and 



/-1.0301 + 0.3060i 
-0.0162 + 0.56181 
0.7134- 0.15041 
y 1.0488 + 0.10861 



-0.3475 - 0.08161 
0.3662-0.14421 
0.2450 - 0.42821 
0.2369 + 0.23461 



According to we obtain that a = 2.01390,6 = 9.84720, c = 0.0566687, d = 0.0125309, r 

(r- 



0.540848, zo = 0.841120. Then, according to Theorem[I]i), we calculate y/l - r"^ zo+( 



-0.1959 < 0, y 1 - r2 + - r) 

Then, we have 



-{cy/bja+d) = 1.8452 > 0. Thus, ^ holds and hence Rs > 0. 



eig(G) = [47.92, -1.46, 0.0016 ± 0.99i, 0.6741, -0.016] . 

There is only one feasible x = 0.6741 among eig(G). Then, we obtain x* = 0.6741 and hence z* 
0.9270. 
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Finally, according to ^ and we have u* = [0.4692 - 0.3024i, 0.1854 - 0.4521i, -0.3258 - 
0.1020i, -0.5655 + 0.1153i]^ and 

Rs = 3.1162 (bits/s/Hz). 

Fig. [T] plots the secrecy rate for different power P. 

Second, we give an example in which Rs > cannot be achieved. We set (e/j,e£;) = (0.05,0.05), 
and hfi = [0.1216 + 0.01181, 0.0106 - 0.0316i, -0.0856 - 0.1063i, 0.2241 - 0.0216i], Iie = [0.3599 + 
0.0174i, 0.1655 -0.19231, -0.2323 -0.40651, 0.7313 -0.22721]. Neither of ^ and ^ holds. Thus, 
Rs > cannot be achieved even when infinite power is used. 

V. Conclusion 

We study the problem of finding the optimal input covariance that maximizes the worst-case secrecy 
rate of a MISO wiretap channel under channel uncertainty, subject to a power constraint. We show that 
the optimal input covariance can be obtained via finding the eigenvalues of a known 6-by-6 matrix. 

Appendix A 
Closed Form Solutions of (fT4l) and (fTsT) 

First, we solve (fTSl) . Let h^; = + eEy- The constraint h^; G Qe becomes ||y|| < 1. Thus, (fT5]) is 
equivalent to 

max ju^h^; + e£;uVP- (31) 

lly||<i 

First, |u1"h£;| + ee is an upper bound of ju^h^; + e^ju^l which can be easily verified from the triangle 
inequality and Cauchy's inequality. Moreover, this upper bound can be achieved: y = (u^h£;/|uth£|)u 
(if u^h^; = 0, then y = u). Thus, the optimal objective value of (OTI ) is (|uth£;| + e^;)^. 

Second , we solve (fT4l) . Let h/j = + e^x. The constraint hji G Or becomes ||x|| < 1. Thus, (fT4l) 
is equivalent to 

mill ju^hi? + e_RU^x|^. (32) 

||x||<l 

Since r* > 1, it holds that the optimal objective value of (l32l) is greater than zero. Note that the 
objective value of ( [32l ) is zero at the point xi = — (u^h/j/e/j)u. Thus, xi must be infeasible and hence 
|u^h/j| > eji- With this, we can show that |u^h/f | —er is a lower bound of |u^h/j + e/fu1"x| which can be 
easily verified from the reverse triangle inequality and Cauchy's inequality. Moreover, this lower bound 
can be achieved: x = — (u^h/j/|u^h/j|)u. Thus, the optimal objective value of (|32l ) is (ju^'h/jj — e/j)^. 
This completes the proof. 
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Appendix B 
Proof of Lemma [U 

First, we show that the optimal objective value is zero if and only if q < 1 — r^. Denote the null space 
of b as E;, with ^^Eb = I and EfeEj^ = I — bb^ Then b^u = if and only if there exists a vector z 
such that u = E;,z which, when inserted into u^aa^^u = q and ||u|p = 1, results in that ||z|p = 1 and 
z^E^aatEfeZ = q. Note that < zTE^aa+EfcZ < Amax(E^aatEb) = atEfoEj^a = at (I - bbt)a = 1 - r^. 
In other words, ||z|p = 1 and ztEj^aa^E^jZ = q hold if and only if g < 1 — r^. The desired result follows. 

Second, we consider the case q > I — r^. The optimal u is a linear combination of a and b, which 
follows from its optimality condition bb^u— /iiaa^u— = or equivalently = [h^ u.)h— {fiia^ u.)a 
where //i and fi2 are multipliers (Obviously, /U2 7^ since r < 1). From this result, we let u = cia + C2b. 
Since e"^u (for any real uj) satisfies the constraints and attains the same objective value as u, we can 
restrict C2 > 0. Inserting u = cia + C2b into the constraints and objective, results in 

|ci|2 + c^ + c5fc2atb + ciC2bta= 1, (33) 
|cip + c^r^ + clc2a^h + cicab^a = q, (34) 
and utbbtu = 1 - |cip(l -r^). (35) 

From (|35] ). we need to maximize |cip. From (l33l ) and (l34b . we get c^l — r^) = I — q which leads to 
C2 = ^/{l — q)/{l — r^) . By denoting ci = |ci|e'^ where 9 is the argument of ci, we can rewrite (1331 ) 
as 

|ci|2 + |ci|c2r(e-^(<^+^) + e^(<^+^)) + (c^ - 1) = (36) 

where (/) is the argument of b^a. This is a quadratic equation of one variable |ci|. It is not difficult to show 
that the optimal 9 = tt — (f>, and the optimal |ci| = C2r + ^/g. Thus, the optimal ci = (c2r + ^)e^^'^~'^\ 
and hence 

Here we have used the fact that, e^^'^~'^^ = —e~^'^ = — a^b/r. Further, from (1351 ). we obtain 

u^bb+u = 1 - (C2r + - r^) 

= (rVg - vT^yr^ )2. (38) 

This completes the proof. 
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P{dB) 

Fig. 1. Worst-case secrecy rate versus power, P, of a MISO wiretap channel with 4 transmit antennas. (e_R, es) ~ (10^^, 10^^). 
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